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■ Abstract 

We calculate (g — 2)/2 of the muon, by improving the determination of the hadronic vacuum 
polarisation contribution, a^ ad,LO , and its uncertainties. The different e + e~ data sets for each 
exclusive (and the inclusive) channel are combined in order to obtain the optimum estimate of 
the cross sections and their uncertainties. QCD sum rules are evaluated in order to resolve an 
apparent discrepancy between the inclusive data and the sum of the exclusive channels. We conclude 
a had,LO _ (5g3 1 -j- 5.9 exp ± 2.0 ra d) x 10 -10 which, when combined with the other contributions to 
(g — 2)/2, is about 3u below the present world average measurement. 



1 Introduction 



The muon anomalous magnetic moment, = (g^ — 2)/2, is one of the most precisely measured 
quantities in contemporary particle physics. The world average of the existing measurements is 

a^ xp = 11659203(8) x 1(T 10 , (1) 

which is dominated by the recent value obtained by the E821 collaboration at BNL[1]. It is so precisely 
measured that it is very useful in probing and constraining New Physics beyond the Standard Model 
(SM). It is therefore important to evaluate the SM prediction of as accurately as possible. 
The SM contribution to a M may be written as the sum of three terms, 

= a« ED + a EW + (2) 

The QED contribution, a^ ED , has been calculated up to and including estimates of the 5-loop contri- 
bution, see reviews [2, 3], 

o qed = n6 584 705.6(2.9) x 1CT 11 . (3) 

In comparison with the experimental error in eq. (1), and with the hadronic contribution error dis- 
cussed later, the uncertainty in a^ ED is much less important than other sources of uncertainty. The 
electroweak contribution a EW is calculated through second order to be [4, 5, 6] 

a EW = 152(1) x HT 11 . (4) 

Here again the error is negligibly small. 

Less accurately known is the hadronic contribution a^ ad . It can be divided into three pieces, 

„had _ „had,LO , ^had.NLO , „had,l-b-l /c\ 
a H — a n ' ix ' fj, • \°) 

The lowest-order (vacuum polarisation) hadronic contribution, a^ ad,LO , has been calculated by a num- 
ber of groups. The value 

a had,LO = 6 9 24( 6 2) x HT 11 (6) 

taken from Ref. [7] has been frequently used in making comparisons with the data. The next-to-leading 
order hadronic contribution, evaluated to be [8, 9] 

a had,NLO = _ 100 ( 6 ) x iQ- 11 . (7) 

The hadronic light-by-light scattering contribution a h&d,\-b-i ^ as ^ een recently reevaluated [10]-[15], 
and it is found to be 

a h ^ d ' l - M = 80(40) x HP 11 , (8) 

where we quote the estimate of the full hadronic light-by-light contributions given in Ref. [16]. From 
eqs. (6), (7) and (8), we can see that aJj ad ' LO has the largest uncertainty, although the uncertainty in 
the light-by-light contribution a h&d,i-M j g a ^ go i ar g e 

In this letter we update the evaluation of a^ ad ' LO , which is given by the dispersion relation 
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where <Tha,d(s) is the total cross section for e + e — > hadrons(+7) at centre-of-mass energy y/s. The 
kernel function -fC(s) is given by 

o N 3S fx 2 , n, (1 +X 2 )(1 +X) 2 / , , X 2 \ 1 + Z n, ] 

> Am D = ^2\y {2 ~ x ) + — ^ — ( ln(1 + x) " x + y ) + t^ x lnx \ ' (10) 

with x = (1 — /3 (U )/(1 + where /3 M = ^/l — 4m^/s. -fT(s > 4m 2 ) increases monotonically from 0.63 
to 1 in the range 4m 2 < s < oo. The form of K for s < 4m 2 is given in [17], and is used to evaluate 
the small ir°j contribution to a^ ad ' L °. 

To evaluate cr had (s), we use experimental data up to 11.09 GeV and perturbative QCD thereafter. 
The most important contribution to o~had(s) comes from the e + e~ — > tt + tt~ channel; the channels 
7T + 7r~7r°, K + K~, KgK®, tt + tt~ tt + ~tt~ , 7r + 7r _ 7r°7r°, etc. give subleading contributions. For example, if 
we evaluate a\ ia d(s) using the sum of the data for the above 6 exclusive channels up to = 1.43GeV, 
we obtain 87% of the total a]f d ' LO , with the vr+vr" contributing 72%. Then if from 1.43 to 11.09 GeV 
we were to use the measurements R(s) = cri m ^(s)/a(e + e~ — > /j + /j~) to determine <Th a d(s), we would 
obtain another 11% of the total. 

As mentioned, our study to update 

a had,LO ig 

motivated by the increasing accuracy of the exper- 
imental value for a^ xp , see eq. (1). The special features of our analysis are (i) that it is data-driven, 
based on all available data, including the new data on exclusive channels from Novosibirsk, partic- 
ularly ir + ir~ [18], and the BES data on the R ratio [19], (ii) the careful application of a clustering 
method, so that data of differing precision can be combined consistently, and (iii) the use of QCD sum 
rules to resolve an apparent discrepancy between the inclusive and exclusive determination of (Jha,d(s) 
for 1.4 < ^ < 2 GeV. 

We chose not to use data on hadronic r decays to further constrain the e + e~ — > 2-7r,4-7r channels 
for y/s < m T , because of the possible uncertainties connected with isospin-breaking effects. A careful 
study of the effects of including r data has been made in Ref . [20] . 



2 Processing the hadronic data 

We apply the hadronic vacuum polarisation corrections given by Swartz [21]. Moreover, we calculated 
the final state radiative effects for all e + e~ — > ir + ir~ data, except for the new CMD-2 data [18], based 
on eq. (45) of Ref. [22]. These two corrections increase the 7r + 7r~ contribution by about 1.1 x 10~ 10 , 
to which we assign a 50% error. For the dominant CMD-2 data the radiative corrections are already 
included by simply taking the cross section numbers for a® w ^y Since the recent accurate data [23] 
in the u and <j> resonance regions have not been corrected for any vacuum polarisation effects (see 
the comment in [20]), we apply the full (including the lepton vacuum polarisation) corrections to 
these data. For the data on the other exclusive channels, and the inclusive data not discussed in [21], 
insufficient information is available to make reliable radiative corrections. We therefore assign an 
additional ±1% uncertainty to their contribution to a^ ad ' LO . The net effect is an error of about 
±2 x 10~ 10 due to radiative corrections; a further discussion will be given in [24]. 

We now come to the important problem of 'clustering' data from different experiments (for the 
same hadronic channel) . To combine all data points for the same channel which fall in suitably chosen 
(narrow) energy bins, we determine the mean R values and their errors for all clusters by minimising 
the non-linear x 2 function 



x\R m , fk) =£[(i- fk) mf + E E [{R\ k ' m} - hRm) /d*f' m} ] • (ii) 

fc=l 771=1 1=1 
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Here R m and ft are the fit parameters for the mean R value of the m cluster and the overall 
normalization factor of the k th experiment, respectively. R\ k ' m ^ and dR\ k ' m ^ are the R values and 
errors from experiment k contributing to cluster m. For di?j fc '"^ the statistical and point-to-point 
systematic errors are added in quadrature, whereas dfk is the overall systematic error of the k th 
experiment. Minimization of (11) with respect to the (-/V e xp + -Wclust) parameters, and R m , gives 
our best estimates for these parameters together with their error correlations. 

Our definition (11) implies piecewise constant R values but imposes no further constraints on the 
form of the hadronic cross section ('minimum bias'). Due to use of the overall normalization factors it 
also results in an adjustment of the different sets within their systematic uncertainties. This means, 
for example, that sparse but precise data will dominate the normalization. Still, the information on 
the shape of R from sets with larger systematic uncertainties is preserved, and all data contribute 
weighted according to their significance. 

The error estimate for each hadronic channel is then done using the complete covariance matrix 
returned by our x 2 minimization. Therefore statistical and systematic (point-to-point as well as 
overall) errors from the different sets are taken into account including correlations between different 
energies (clusters) . The minimum x 2 directly reflects the quality of the fit and the consistency of the 
data. We have checked that for all hadronic channels we find a stable value and error for a^ ad,LO , 
together with a good 1 x 2 if we var T the minimal cluster size around our chosen default values (which 
are typically about 0.2 MeV for a narrow resonance and about 30 MeV for the continuum). 

The dispersion integral (9) is performed integrating (using the trapezoidal rule) over the clustered 
data directly for all hadronic channels, including the to and (j) resonances. Thus we avoid possible 
problems due to missing or double-counting of non-resonant backgrounds, and interference effects are 
taken into account automatically. As an example we display in Fig. 1 the most important ir + ir~ 
channel, together with an enlargement of the region of p-u interference. 

In the region between 1.43 and ~ 2 GeV we have the choice between summing up the exclusive 
channels or relying on the inclusive measurements from the 772, MEA, M3N and ADONE experiments 
[25] . Surprisingly, the sum of the exclusive channels overshoots the inclusive data, even after having 
corrected the latter for missing two-body and (some) purely neutral modes. The discrepancy is shown 
in Fig. 2, where we display data points with errors after application of our clustering algorithm. 



3 Results 

In Table 1 we list the contributions to a^ ad,LO from different energy regimes: From the two-pion 
threshold up to 0.32 GeV, chiral perturbation theory is applied (see e.g. [26]-[28]) and for the high 
energy tail above 11.09 GeV, R is calculated using perturbative QCD. The J/ip, ip(2S) and Y(15 — 65) 
resonance contributions are evaluated in narrow-width-approximation. Apart from those contributions 
we use the direct integration of clustered data as described above. For the controversial region from 
1.43 to 2 GeV we present two results: if we use the lower lying inclusive data the corresponding 
contribution is considerably smaller than the one resulting from the use of the sum over the exclusive 
channels. 

A more detailed breakdown of the contributions, and a full discussion of the data used, will be 
presented elsewhere [24]. There, we will also present an updated value of the QED coupling at the Z 
pole, a(M|). 

1 In the channels e + e~ — > tt + ir~ tt + ir~ , tt + -k~-k°, in which data sets are mutually incompatible, Y 2 /degree of freedom 
= 2.1, 1.6. For both cases the error is enlarged by a factor of \J x 2 /dof . 



3 



3 1400 
O 

1200 
1000 
800 
600 
400 
200 




0.3 




0.4 0.5 0.6 0.7 0.8 0.9 



1.1 1.2 
Vs (GeV) 



^ 1500 
^ 1400 - 

1300 

1200 

1100 - 

1000 ^ 
900 r 
800 E" 
700 E~ 
600 

0.74 



0.75 




0.76 0.77 0.78 0.79 



0.8 



0.81 0.82 
Vs (GeV) 



Figure 1: e + e — > 7r + 7r data up to 1.2 GeV, where the shaded band shows the result of clustering 
The second plot is an enlargement of the p-to interference region. 
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Figure 2: The inclusive and the sum of exclusive channel values of R, after the data have been 
clustered. 



4 Resolution of the ambiguity: QCD sum rules 

To resolve the ambiguity between the exclusive and inclusive data values for R(s) in the range 1.43 < 
y/s < 2 GeV (see Fig. 2), we evaluate QCD sum rules of the form 

f ° ds R(s)f(s) = [ ds D(s)g(s) (12) 

Js th JC 

where sq is chosen just below the open charm threshold and C is a circular contour of radius sq. D(s) 
is the Adler D function, 



ds \ s 



where R(s) 



ImII(s). 

s 



(13) 



We use the experimental data 2 for R(s) (or equivalently crhad(s)) to evaluate the left-hand-side, while 
QCD is used to determine [29] 

D(s) = D (s) + D m (s) + Aip(s), (14) 

where Dq is the 0(a|) massless, three-flavour QCD prediction, D m is the (small) quark mass correction 
and D np is the (very small) contribution of the condensates. We take f(s) to be of the form (1 — 
s/so) n (s/so) m , with n + m = 0, 1 or 2. Once f(s) is chosen, the functional form of g(s) is readily 
evaluated. For example, the n = 1, m = sum rule is 



ds 



(15) 

/s t h "0/ 27T V 2s 2s/ 

The sum rules with n = 1 or 2 and m = are found to maximize the fractional contribution of the 

left-hand-side of (12) coming from the relevant 1.43 < yfs < 2 GeV interval. The evaluations of 

these two sum rules are shown in Table 2. Consistency clearly selects the inclusive, as opposed to 

the exclusive, determination of R(s). A more detailed discussion of the QCD sum rules, and their 

2 The J/ip and tp(2S) resonance contributions are, of course, omitted. 
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energy ran 


?e (GeV) 


comments 


a had,LO x 1Q 10 

f- 1 


2m^ . . 


0.32 


chiral PT 


2.30 ±0.05 


0.32.. 


1.43 


excl. only 


596.73 ±5.18 


1.43.. 


2.00 


excl. only 


38.14 ± 1.72 






incl. only 


32.43 ± 2.46 


2.00... 


11.09 


incl. only 


42.09 ± 1.25 


J/Y> and V(25) 


nar. width 


7.31 ± 0.43 


T(15- 


-65) 


nar. width 


0.10 ±0.00 


11.09. 


. . oo 


pQCD 


2.14 ±0.01 


E of all 


ex-ex-in 
ex-in-in 


688.81 ±6.17 
683.11 ±5.89 



Table 1: A breakdown of the contributions to different intervals of the integration (9) for ' . The 
alternative numbers for the interval 1.43 < ^fs < 2 GeV correspond to using data for either the sum 
of the exclusive channels or the inclusive measurements, see Fig. 2. The total also includes a small 
0.13 x 10~ 10 contribution from the 7r°7 channel near its threshold (also included in the second line 
above). 



sum rule 


l.h.s. (data) 


r.h.s. (QCD) 


n = 1, m = 


15.34 ± 0.39 (incl) 
15.99 ± 0.35 (excl) 


15.34 ±0.08 


n = 2, m = 


10.40 ± 0.25 (incl) 
10.90 ± 0.22 (excl) 


10.30 ±0.06 



Table 2: The results of evaluating sum rule (15) and the corresponding one with f(s) = (1 — s/sq) 2 , 
where ^ = 3.7 GeV. The main QCD error comes from as(Mf) = 0.117 ±0.002 [30]. The 'incl' and 
'excl' alternatives refer to using the inclusive or exclusive e + e~ data in the region 1.43 < \/s < 2 GeV, 
see Fig. 2. 

evaluation, will be given in [24]. 

The same conclusion with regard to the resolution of the inclusive/exclusive ambiguity in the 
1.43 < t/s < 2 GeV interval was reached in an independent analysis [31]. 

5 Conclusions 

We have undertaken a data-driven determination of the hadronic vacuum polarisation contribution to 
a had,LO_ yy e j^yg usec j a ll available e + e~ data and a non- linear \ 2 approach to cluster data for the 
same channel in narrow bins. In particular, the fit allows the normalizations of the individual data 
sets to be collectively optimized within their uncertainties. We found that there was a discrepancy 
between the inclusive value for a(e + e~ — ► hadrons) and the sum of the exclusive channels in the region 
1.4 < yfs < 2 GeV, which gave an uncertainty of about 6 x 10~ 10 in a^ ad ' L °. We used a QCD sum 
rule analysis to resolve the discrepancy in favour of the inclusive data. Thus finally we find that the 
SM predicts 

had.LO = ( 68 3.i ±5 .9 cxp±2 .0 rad ) x 10- 10 . (16) 
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Summing up all SM contributions to a^ M as given in eqs. (2)-(8), with (6) replaced by (16), we 
conclude that 

a™ = (11659166.9 ± 7.4) x 10~ 10 , (17) 

which is 36.1 x 1CP 10 (3.3a) below the world average experimental measurement. If, on the other 
hand, we were to take the value of a^ ad ' LO obtained using the sum of the exclusive data in the interval 
1.43 <^s<2 GeV then we would find a™ = (11659172.6 ±7.7) x 1(T 10 , which is 30.4 x 1(T 10 (2.7a) 
below af, xp . 

An independent SM prediction has very recently been made [20]. Their final e + e~-based result, 
(684.7 ± 6.0 CX p i 3.6 ra d) x 10~ 10 , is very similar to ours. However, the overall agreement hides larger 
differences in individual contributions (but within the quoted uncertainties). Our result (16) agrees 
also fairly well with a recent reevaluation of the leading hadronic contribution from F. Jegerlehner, 
who also used the recent CMD-2 data [18] and obtained (688.9 ± 5.8) x 10~ 10 , see [32]. In order to 
facilitate a comparison with these two predictions, we will present a detailed breakdown of our result 
elsewhere [24]. 

For the future, we can expect further improvement in the accuracy of the experimental (g — 2)/2 
measurement. As far as the SM prediction is concerned, we may anticipate low energy data for a 
variety of e + e~ channels, produced via initial state radiation, at the ^-factory DA<I>NE [33] and at the 
£?-factories, BaBar and Belle, see, for example, [34]. For instance, by detecting the it + it~j channel, 
it may be possible to measure e + e~ — ► ir + ir~ as low as about yfs = 400 MeV. Moreover, CMD-2 
measurements of e + e~ — > tt + tt~ have already been made in this low energy region [35]. When these 
latter data are final, we anticipate that they would improve the error on aJj ad ' LO by about 1 x 10~ 10 . 
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